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SUMMARY 

Heat-transfer rates, velocity profiles, and temperature profiles 
for the turbulent incompressible flow of air over a flat plate with a 
constant surface temperature have been measured at Reynolds numbers up 

to 3.5X10 • The turbulent heat- transfer measurements agree well with 
the von Karman analogy, and the velocity profiles agree with the data of 
previous investigators. The temperature profiles are similar to the ve- 
locity profiles, both being adequately described by power formulas. 

INTRODUCTION 


The present report is the first of a series of four covering a 
three-year investigation of heat transfer in. the turbulent incompressible 
boundary layer with arbitrary surface temperature (see ref. l) . In this 
report the experimental apparatus is described and the results of ex- 
periments with constant surface temperature are presented. Part II con- 
tains experimental results and analyses for a step' temperature distri- 
bution (ref.' 2). In part III the step-function analysis is used to pre- 
dict heat-transfer rates for several variable -wall- temp’erature cases, and 
the predictions are compared with experiment (ref. 3) . A simple method 
for handling variable-surface-temperature problems is presented. The 
effect of the location of transition on the heat transfer in the turbu- 
lent boundary layer is analyzed and compared with experiments in part IV 
(ref. 4) . 

The broad objectives of this program were to investigate experimen- 
tally the problems of heat transfer in a turbulent incompressible bound- 
ary layer on a flat plate, with negligible pressure gradient, at high 
Reynolds numbers. The Mach number, and temperature difference are suffi- 
ciently low that compressibility effects are negligible, and temperature- 
dependent fluid-property effects are small. This is a problem that has 
been extensively investigated analytically, but adequate experimental con- 
firmation of the analyses was lacking. 
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This first report treats the problem of heat transfer from a flat 
plate at constant wall temperature. A number of analyses for the heat- 
transfer rate have been proposed, and some experiments have been per- 
formed. However, the only reliable data at Reynolds numbers over about 
5 

5X10 have been obtained at high velocities, where compressibility ef- 
fects are important. Thus, there is a need for experimental confirma- 
tion of the analyses, and the constant-surface-temperature data con- 
tained in this report are presented primarily for purposes of this con- 
firmation. They also form a point of departure for the work described 
in the subsequent reports. 

Perhaps the best analysis available at the present time is the heat- 
transfer - momentum analogy of von Karman (ref. 5) , by which the local 
heat-transfer coefficient has been related to the local friction factor. 

A suitable expression for the friction factor in terms of the local Reyn- 
olds number is then required. A number of analyses for friction factor 
are available, and that of Schultz-Grunow (ref. 6) has been experimen- 
tally verified over a wide range of Reynolds numbers. The von Karman 
analogy, combined with the friction-factor formula of Schultz-Grunow, 
allows prediction of the local heat transfer. This result is in excel- 
lent agreement with the experimental data to be- presented in this report 

5 S 

at flow Reynolds numbers from 10 to 3.5X10 . 

An alternative method of analyzing skin friction and heat transfer 
is to assume the form of the velocity and temperature profiles, and then 
to use the momentum and energy integral equations of the boundary layer 
to arrive at expressions for the friction factor and heat-transfer coef- 
ficient as functions of the local Reynolds number. In making such anal- 
yses it is convenient to represent the profiles by power expressions, and 
it is commonly assumed that velocity and temperature vary as the one- 
seventh power of the distance from the surface. The data of Schultz- 
Grunow indicate that the l/7-power velocity profile is reasonable at 

7 

Reynolds numbers near 10 but is not accurate for lower or higher Reyn- 
olds numbers (ref. 6) . Velocity profiles obtained in the present inves- 
tigation agree quite well with the Schultz-Grunow profiles; the temper- 
ature profiles obtained at constant wall temperature are indeed "similar" 
to the velocity profiles, both being described very well by l/5. 6-power 
formulas. The survey data presented in this report are therefore of in- 
terest because they may be used as the basis for integral analyses for 
friction and heat transfer in the turbulent boundary layer. 

This investigation was carried out at Stanford University under the 
sponsorship and with the financial assistance of the National Advisory 
Committee for Aeronautics. The assistance of the following people is 
gratefully acknowledged: H. M. Satterlee designed the experimental ap- 

paratus and supervised the construction with B. J. Grotz, who constructed 
the instrumentation and assisted in some of the preliminary tests. 
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P. W. Rundstadler, R. M. Foster, J, D. Stephenson, G. L. Meyers, and 
M. L. van der Ploeg assisted in the wind-tunnel testing, and L. R. 
Reneau and H. Singh reduced the hulk of the data. 



G 

h 

k 

Pr 

P 

P S 


Pt 

q hl 

Re x 

Re 5 


St 


At 


SYMBOLS 

friction factor, T w / (pu£/ 2) 
pressure coefficient, (p - P s )/ (p-t “ P s ) 
specific heat at constant pressure, Btu/(lb)(°F) 
free- stream mass velocity, pu , lb/(hr) (sq ft) 

convective heat-transfer coefficient, q^./At, Btu/(hr) (sq ft) (°F) 
thermal conductivity of fluid, Btu/ (hr) (ft) (°F) 

Prandtl number, pc /k 

tT 

pressure, lb/sq ft 

static pressure upstream of plate > lb/sq ft 
total pressure upstream of plate, lb/sq. ft 
heat flux in boundary layer, Btu/ (hr) (sq ft) 
heat flux at wall, Btu/ (hr) (sq ft) 
flow Reynolds number, Gx/p 
Reynolds number based on 8, GS/p 
local Stanton number, h/GCp 
absolute wall temperature, °R 
absolute free-stream temperature, °R 

k - °t 

temperature in boundary layer, °F 
mean temperature of heated strip, °F 
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€ H 

€ m 
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M* 
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p 


- V ° F 

wall temperature , °F 
free -stream temperature, °F 

dimensionless temperature, (t w - t^ ) pCp-y/ x w /p/ q^ 

velocity in x-direction, ft/ sec 

free-stream velocity, ft/sec 

dimensionless velocity, u/ ->J T w /p 

velocity in y-direction, ft/sec 

distance from leading edge, ft 

distance from plate, ft 

dimensionless distance from wall, x w /p/v 
thermal diffusivity of fluid, sq ft/hr 
thickness of hydrodynamic boundary layer, ft 

displacement thickness, n - < 
thickness of thermal boundary layer, ft 


«/<J 


dy, ft 


conduction thickness, 



0)dy, ft 


eddy diffusivity for heat, sq ft/hr 
eddy diffusivity for momentum, sq ft/hr 
dimensionless temperature, (t w - t-^ ) / (t w - t^,) 
viscosity of fluid, lb/ (hr) (ft) 
kinematic viscosity, |i/p, sq ft/hr 
fluid density, lb/cu ft 
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shear stress in boundary layer , lb/sq ft 
t w shear stress at wall, lb/ sq ft 

ANALYSIS 
Heat Transfer 


The problem of turbulent heat transfer from a flat plate has been 
attacked analytically by a number of investigators, generally by one of 
two methods. One approach is the integral method, wherein the form of 
the velocity and temperature profiles is assumed, and the energy inte- 
gral equation is used to arrive at a relation between the local heat- 
transfer coefficient and the local friction factor. The integral method 
is very powerful, because the result is relatively insensitive to the 
choice of the temperature and velocity profiles. The second method is 
the analogy method, wherein mechanisms for heat and momentum diffusion 
throughout the boundary layer are assumed, and an empirical velocity 
profile is employed to determine the heat-transfer rate in terms of the 
local friction factor. The resulting expressions for the local heat- 
transfer coefficient are generally more complicated than those obtained 
by integral methods, but they agree better with experimental values. 

It appears that the best analysis available at the present time is 
the heat-transfer - momentum analogy of von Karman (ref. 5), who 
obtained 



For Prandtl numbers of unity this result reduces to the familiar Reynolds 
analogy, 


St = C f /2 (2) 

It is evident that the friction factor must be known if the heat- 
transfer rates are to be determined, and a number of friction- factor 
analyses have been made (ref. 7). An expression for the local friction 
factor may be derived by integral methods, if it is assumed that the ve- 
locity profile follows a l/7-power law and that the wall shear stress is 
related to the boundary-layer thickness by 

Cj/2 = 0.0228 Reg 1 / 4 


(3) 
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This relation was first proposed by Blasius on the basis of pipe-flow 
friction data. These assumptions lead to (ref. 8) 

C f /2 = 0.0296 Re^ 0, 2 (4) 

Equation (4) appears to be adequate at Reynolds numbers below 10® but 

r 7 

gives friction factors that are low above Re v of 10° to 10 . A more 

refined analysis was made by Schultz-Grunow, who measured velocity pro- 
files and local friction factors at Reynolds numbers up to 10 9 (ref. 6) . 
He then used the momentum integral equation to find the friction factor 
and obtained a result that may be represented by 

C f /2 = 1. 60 (in Re x )" 2 * 58 (5) 

Equation (5) is in good agreement with Schultz-Grunow ' s data for Reynolds 

Q 

numbers up to 10 , and it is felt that this is the best friction- factor 
relation available at the present time. 

The von Karman analogy (eq. (l) ) and the Schultz-Grunow friction 
formula (eq. (5) ) may be combined to give the local heat-transfer coef- 
ficient in terms of the local Reynolds number: 

1. 60(ln ReJ" 2 * 58 

St = — " ~ ~r =r (6) 

1 + 1.2 6 (in. Re x ) _1 * 29 [5Pr + 5 ln(5Pr + l) - 14J 

For air, which has a Prandtl number -around 0.7 (see ref. 7), equation , 

(6) may be represented approximately in the range 10® < Re x < 10^ by 

StPr 0,4 = 0.0296 Re" 0,2 (7) 

This relation is easier to use for calculations than equation (6). Com- 
parison of equations (7) and (4) shows that, in this range, 

StPr 0,4 = C f /2 (8) 

Equation (8) represents a modification of the familiar Colburn analogy 
(ref. 9) : 

StPr 2 / 3 = C f /2 

The Colburn relation predicts heat- transfer coefficients for air that are 
too high, and the modification (eq. (8) ) is better. 
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The foregoing equations were obtained only for constant fluid prop- 
erties. If the temperature of the plate is considerably different from 
that of the free stream, there may be considerable variation in the fluid 
properties appearing in the Stanton and Reynolds numbers, and a question 
arises as to the best temperature for evaluation of these properties. 

It can be shown (ref. 10 ) that an adequate method for taking into consid- 
eration the influence of temperature -dependent fluid properties for gas 
flow in both internal and external boundary layers is to evaluate all 
properties at the free-stream static temperature and then to include all 

properties in a factor (T^/T^) 131 , where the exponent m is a function of 

geometry alone. Examination of the results of reference 11 indicates 
that, for the turbulent incompressible boundary layer, the Stanton number 

varies as (T^/T^) other things being equal. This observation may be 

used to correct the foregoing equations for temperature -dependent fluid- 
property effects. Thus, for the simple power relation (eq. (7)) one may 
write 


-0.4 


StPr 0 ’ 4 = 0.0296 ReJ 0 * 2 ^ 



where the Stanton number and Reynolds number are to be evaluated at the 
free-stream temperature . 


Velocity Profiles 


Data correlation . - Velocity survey data in turbulent boundary lay- 
ers are usually correlated in one of three ways: , 


u /h» = fi(y/s) 

u + = f 2 (y + ) 

v - u+ = f 3 Cy/ 8 ) 


( 10 ) 

( 11 ) 

( 12 ) 


where u + is a dimensionless velocity, defined as 


u + - 


u 




and y is a dimensionless distance, defined as 

yV T w/p 


y + = 


V 
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A plot of the type (10) may be prepared simply from measurements of 
the velocity in the boundary layer as a function of the distance from the 
wall. The methods (ll) and (l2) require additionally the determination 
of the wall shear stress t . Method (lO) is most satisfactory in the 

outer portions of the boundary layer, where the flow can be characterized 
by the thickness of the hydrodynamic boundary layer 5. Von Karman .(ref. 
12) has given theoretical basis for the "universal velocity profile," 
used in the method of (ll) . However, it has been found experimentally 
that this method of correlating profiles works best near the wall, where 
free- stream conditions have the least influence. The third type of plot 
(12) , which is often referred to as the "universal velocity deficiency 
law," represents attempts to tie the wall effects to the free-stream ef- 
fects by introducing both the wall shear stress and the boundary-layer 
thickness. Data plotted on this basis correlate nicely in the outer re- 
gions of the boundary layer (y/& > 0.01) , but the method fails near the 
wall. 


The boundary-layer "thickness" is a rather nebulous thing, since in 
reality the boundary-layer velocity reaches the free-stream velocity at 
an infinite distance from the plate. However, the velocity is almost 
equal to the free-stream velocity a short distance from the plate ; and 
it is this distance that is usually referred to as the "boundary-layer 
thickness." There then arises a question as to how much is "almost," or, 
in other words, what fraction of the free-stream velocity occurs at the 
"edge" of the boundary layer. Often the distance at which the velocity 
is 99 percent of the free-stream velocity is taken as the boundary-layer 
thickness; however, this definition has no physical meaning and moreover 
is difficult to determine accurately from experimental data. In integral 
treatments of the boundary layer, some relation of the type (10) would be 
advantageous, where a "boundary-layer thickness" is used to characterize 
velocity profile in the boundary layer. Therefore, some means is desired 
of evaluating for experimental surveys the same boundary-layer thickness 
used in the integral methods. A meaningful boundary-layer thickness can 
be determined by making use of the fact that turbulent velocity profiles 
can, to a good approximation, be represented by equations of the form 



The power parameter m is about 5 to 8. For a profile of this type, the 
boundary-layer thickness 5 is related to the displacement thickness 
5 , which has real physical meaning, by 

5 - (l + m)S* (14) 

The displacement thickness 5* can be determined quite accurately by 
integration of the velocity profile; m can be determined by plotting 
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u/uo, against y on log-log paper. Then, 8 may he calculated from 

equation (l4) . This technique allows calculation of a meaningful 8 
that corresponds to the 8 used in integral analyses. This method has 
been used in the reduction of velocity survey data obtained in this in- 
vestigation; the values of 8 calculated in this manner are quite close 
to the values of the "99-percent" boundary-layer thickness. 

Universal velocity profile for flat plate . - A velocity profile of 
the type (11) is well known for fully established turbulent flow in a 
pipe; this is the familiar universal velocity profile for turbulent flow 
in pipes, which is based largely on the experimental data of Nikuradse 
(ref. 8). The velocity profile in an external turbulent boundary layer 
is not expected to be too different from that found in a pipe, and this 
assumption has been used by numerous investigators in analyzing turbu- 
lent boundary layers (ref . 5) . The main differences between internal and 
external shear flows occur away from the wall, where the free-stream 
conditions may influence the boundary layer. However, the two flows 
should be quite similar near the wall, where the wall effects predominate. 
A great deal of accurate flat-plate velocity-profile data has been ob- 
tained by Schultz -Grunow (ref. 6), who also made wall shear-stress meas- 
urements. However, these data are limited to the outer regions of the 
turbulent boundary layer, and no data were obtained in the laminar sub- 
layer region. However, by combination of the Schultz -Grunow flat-plate 
survey data and the universal velocity profile for pipe flow, a suitable 
universal velocity profile can be constructed for a flat plate. This has 
been done here, and the details of this combination are presented in the 
following paragraphs. 

The boundary layer has been divided into four distinct regions : 

(1) A laminar sublayer formed near the wall 

(2) A buffer region adjacent to the laminar sublayer 

(3) A turbulent core adjacent to the buffer layer and extending 

over about l/3 of the boundary layer 

(4) A turbulent wake, which extends from the core to the free stream, 

approximately 2/3 of the boundary layer 

In the laminar sublayer, the mechanism for momentum diffusion is 
entirely viscous, so that 

du 

T bZ = ^ 

Because this sublayer is so thin, the shear stress is essentially con- 
stant and equal to its value at the wall Thus, by suitable 



10 


manipulations and integration, the preceding relation leads to the equa- 
tion for the velocity in the sublayer: 

u + = y + (15) 

Very little data have been obtained in the laminar sublayer; therefore , 
it is difficult to tell how far from the wall it extends. It is common 
to assume that, in a pipe, the sublayer extends to y + = 5, and it is 
assumed that this value is reasonable for the external boundary layer as 
well. 

The buffer layer is the region near the wall where momentum diffu- 
sion by turbulent eddies becomes more important. Velocity surveys in 
pipes indicate that the buffer layer may be described by (ref. 5) 

u + = -3.05 + 5.0 In y + (16) 

The extent of the buffer layer can be examined by plotting u + against 
y + on semilogarithmic paper, which renders (16) a straight line (see 
fig. l) . At the outer edge of the buffer layer, a distinct break is ob- 
served, and in the pipe-flow data this break occurs at about y + = 30. 
Upon examination of these data, von Karman took y + = 30 as the outer 
edge of the pipe -flow buffer layer (ref. 5) . The data of Schultz -Grunow 
do not extend below y + = 50 (ref. 6), but an extrapolation of his data 
towards the wall intersects the pipe-flow buffer layer (eq. (16) ) at 
about y + = 18.2. It therefore seems reasonable to give this value as the 
outer edge of the buffer layer for the flat-plate turbulent boundary 
layer. 

The velocity profile in the turbulent core of the flat-plate bound- 
ary layer is obtained entirely from Schultz -Grunow' s data, which indi- 
cate that, in the core, 

u + = 4.4 + 2.43 In y + (17) 

Figure 1 indicates that the core extends from -y + of 18.2 to 360. 

In the region of the turbulent wake, for y + greater than 360, 
Schultz-Grunow’s data exhibit considerable scatter, and the best inter- 
pretation of the data in this region appears to be 

u + = -4.4 + 3.96 In y + -' (18) 

Because of the influence of free-stream conditions, the velocity data do 
not correlate as well in the turbulent wake as in the inner regions, as 
both the data of Schultz -Grunow and the present data show (fig. l). 
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To summarize, a universal velocity profile for the flat plate has 
been obtained by "patching" the boundary-layer data of Schultz -Grunow to 
the familiar universal velocity profile found for turbulent flow in pipes 
This results in a four-region boundary layer as follows: 

(1) Laminar sublayer (eq. (15 ) ) : 0 < y + < 5 

(2) Buffer layer (eq. (16) ) : 5 < y + < 18,2 

(3) Turbulent core (eq. (17) ) : 18. 2 < y + <: 360 

( 4 ) Turbulent wake (eq. (18)) : y + > 360 


Temperature Profiles 


Similarity of velocity and temperature fields . - Because of the 
similarity between the mechanisms of heat and momentum transfer, the tem- 
perature field in a boundary layer can be examined by consideration of 
the velocity profile. This is especially true for a flat plate at con- 
stant temperature, as can be seen by examination of the differential 
equations and boundary conditions for heat and momentum transfer. The 
momentum equation for the turbulent incompressible boundary layer on a 
flat plate may be written as 


u 

u 


cbc 




d(u/Uo) 

5y~ 


— ? 

^ ~5y 


(v + e M ) 


<^y 


(l9a) 


where the shear stress in the boundary layer is given by 


M = (v + e M ) M % 2 ll (19b) 

P ay 

Equation (l9a) is subject to the boundary conditions 

n(x,0) _ Q u(x,q°) _ ' 

Udo 


Similarly, the energy equation of the turbulent incompressible boundary 
layer may be written as (for constant wall temperature, dissipation terms 
neglected) 


u 50 v 50_ 1 5 
u^ + u^-u^ 



(20a) 


where the heat flux in the boundary layer is 
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^bl / . bl 

~ - (a + E H> UjT 

The dimensionless temperature 0 is defined by 

** " Hi 


9 = 


t - t 

V GO 


= a(x,y) 


(20h) 


The energy equation is subject to the boundary conditions 

©(x,0) = 0 0(x,oo) =s 1 

The quantities and eg are defined as the eddy diffusivities 

for momentum and heat, respectively, and it is well known that the ratio 
€ h/ € M -^ s near unity for gases and other fluids with Prandtl numbers near 

1. If the Prandtl number is 1, and if the eddy diffusivities are equal, 
the momentum and energy equations of the boundary layer are identical in 
form. Moreover, since the plate is at constant temperature, their boundary 
conditions are identical. Therefore, the solutions of the two equations 
must be the same, and the velocity profile u/u^ is equal to the 
dimensionless-temperature profile 6. 

If the fluid has a Prandtl number near unity, and the diffusivities 
are approximately equal, the temperature and velocity profiles would be 
expected to be quite similar in shape. It is generally observed that, 
for the turbulent boundary layer of air on a flat plate at constant tem- 
perature, the velocity and temperature profiles are similar when based 
on their own boundary-layer thicknesses. In other words, experiments 
indicate that , 



where S and & T are the thicknesses of the hydro dynamic and thermal 

boundary layers, respectively. This similarity is of great importance 
in integral treatments of turbulent heat transfer. Rube sin (ref. 13) 
shows that, for constant wall temperature, the hydrodynamic and thermal 
boundary layers are related by 

h = Pr- ? /l2 

Thus, for air, having a Prandtl number of 0.7, the thermal boundary layer 
is about 23 percent thicker than the hydrodynamic boundary layer. It is 
important to emphasize that similarity is obtained only when the profiles 
are characterized by their own boundary-.layer thicknesses. 
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Universal temperature profile . - An alternative method of deter- 
mining the temperature profile in the turbulent boundary layer for fluids 
with Prandtl numbers different from unity is to use the analogy between 
heat and momentum transfer. This idea was first suggested by von Karman 
(ref. 5), A dimensionless temperature t + may be defined as 

t+ = ^w - 


By dimensional analysis, one may reason that t + , like u + , should be a 
function of y + only (see ref. 12) . Under this assumption, equations 
(19b) and (20b), respectively, may be written as 


and 



If one has on hand a suitable universal velocity profile (u + against 
y + ), postulates the Reynolds analogy (eg * e^) , and makes some assump- 
tions about the shear-stress and heat-flux distribution in the boundary 
layer, he may compute from equations (21) and (.22) the universal temper- 
ature profile, t + against y + . This will now be done for the flat- 
plate profile discussed previously. 


In the laminar sublayer, heat and momentum transfers are due en- 
tirely to molecular transports, and turbulent transfers are unimportant. 
Thus, the eddy diffusivities for heat and momentum are both zero. More- 
over, the sublayer is so thin that the shear stress and heat flux are 
essentially constant through the layer. Thus, by dividing equation (22) 
by (21) and noting from (15) that u + — y + , the following is obtained: 

dt + = Pr dy + 

which may be integrated to give the temperature distribution in the sub- 
layer (since t + = 0 when y + = 0) , 

t + = Pr y + 0 < y + < 5 (23) 

In the buffer layer, the shear stress and heat flux are again as- 
sumed to be constant. However, turbulent heat and momentum transfer be- 
gins to be important in this layer, so that both molecular and turbulent 
effects must be considered. The eddy diffusivity for momentum may be 
determined by substituting the velocity profile (eq. (16) ) into (21) and 
integrating from y + = 5, where u + = 5: 
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Now postulating the Reynolds analogy/ that e H = e M , equation (22) may 
he integrated from y + = 5, where t + = 5 Pr, to obtain 

t + = 5Pr + 5 ln^Pr ^ + 1 - Pr^ 5 < y + < 18.2 (25) 

In the turbulent core the laminar terms are negligible, as practi- 
cally all heat and momentum are transferred by turbulent eddies. If it 
is assumed that the shear stress and heat flux vary in the same manner 
in the outer regions of the boundary layer, then division of equation 
(22) by equation (21) gives 

dt + = du + (26) 

Integrating this relation from y + = 18.2, where t + = 5Pr + 5 ln(2.65Pr 
+ l) , and u + = 11.45, and substiting the velocity profile (17) result 
in 

t + = 5Pr - 7.05 + 5 ln(2.64 Pr + l) + 2.43 In y + 18.2 < y + < 360 (27) 

In the turbulent wake, the laminar terms are again neglected, and 
it is assumed that the shear stress and heat flux vary in the same man- 
ner. Then, using equation (26) and the velocity profile equation (18) , 

t+ = 5Pr - 16.1 + 5 ln(2.64 Pr + l) + 3.96 In y + 360 < y + < 6* (28) 

Naturally, equation (28) is valid only inside the thermal boundary layer, 
where y + is less than 5* = 8 t V\7p/ v ’ 

It should be emphasized that the -universal temperature profile de- 
rived applies only to a flat plate at constant temperature. Only in this 
case are the differential equations and boundary conditions for heat and 
momentum transfer similar in form, which makes the assumption that the 
shear stress and heat flux vary in the same manner reasonable. Figure 2 
shows the predicted temperature profile compared with some experiments 
of the present investigation. 


EXPERIMENTAL APPARATUS AND PROCEDURE 
Apparatus 

The plate used in this investigation had an active flow length of 
60.5 inches and was tested in the 7. 5-foot -diameter free- jet wind tunnel 
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at the Guggenheim Aeronautical Laboratory of Stanford University. Reyn- 
olds numbers up to 3.5X10° could be obtained with air velocities up to 
130 feet per second. 

The active surface was built up of 24 individually heated copper 
strips , which were thermally -insulated from each other. The surface was 
varnished and rubbed down several times to give a hydraulically smooth 
surface. The inactive side of the plate was well insulated so that the 
"back leak" was minimized. Heat meters were installed at several places 
on both the inactive side and at the ends of the heated strips so that 
the back leak and "end leak” could be measured accurately. The emissivity 
of the heated surface was measured in order that radiation from this sur- 
face could be determined. The insulation between the various strips pro- 
vided a "heat meter/' so that conduction between strips could be 
estimated. 

Iron-constantan thermocouples were located near the surface at the 
center of each strip and at several other points in the strip. The ther- 
mocouples were referenced to a small copper plate mounted on one side of 
the heated surface; the reference plate in turn was referenced to 
distilled-water ice at 32° F, This arrangement provided for a direct 
measurement of the temperature difference. Pressure taps located at 
several points on the active surface allowed measurement of static pres- 
sure on the flat plate. A wattmeter was used to determine the power dis- 
sipated by each heater. The plate construction is described in detail 
in reference 14 and is considered at greater length in appendix A of 
reference 1. 

Because of the end leak and because the heaters did not extend the 
full width of the strips, the center temperature was somewhat higher than 
the average temperature. The mean temperature was analyzed in terms of 
the center temperature, which allows correction of the center temperatures 
to obtain the correct mean temperature differences. The analysis was 
verified by measurement of the transverse temperature distribution, which 
was then integrated to give the mean temperature; good agreement with the 
analysis was obtained, the correction being about 13 percent. The anal- 
ysis of transverse temperature distribution is described in appendix B 
of reference 1. 

The plate is shown before installation in the wind tunnel in figure 
3(a). The dimensions of the active side are given in figure 3(b). The 
heated strips are wider nearer the front in order to minimize contamina- 
tion of the hot boundary layer by sidewise mixing with the cold boundary 
layers on the unheated portions. 

Figures 4(a) and (b) show the plate installed in the wind tunnel. 

The thermocouple, power, and pressure leads can be seen coming from the 
sides of the plate. In order to facilitate handling of probes near the 
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active surface , a traversing mechanism was constructed. This mechanism 
(fig. 4(c)) consisted of a streamlined beam suspended on tracks beneath 
the plate. The probes were rigidly supported by a micrometer arrange- 
ment well in front of the beam and could be moved normal to the surface 
or in the flow direction. 

The free- stream velocity was measured with a Pitot-static tube. Ve- 
locity surveys in the boundary layer were made with a specially con- 
structed Pitot tube made of a hyper demnic tube having an outside diameter 
of 0.020 inch and an inner diameter of 0.010 inch. Temperature surveys 
were made with a specially constructed thermocouple consisting of a butt- 
welded iron-constantan thermocouple flattened to a thickness of 0.002 
inch and supported by heavier iron and constantan wires. The thermo- 
couple probe was referenced directly to the free-stream temperature. 

This probe is shown in figure 5. 


Data-Recording Procedure 

To make a test run, the power to each strip was adjusted to give the 
desired wall -temperature distribution* After steady-state conditions 
were obtained, the center temperature and the power to each strip were 
measured. The free-stream velocity over each strip was then measured 
with a Pitot-static tube that was positioned by the traversing mechanism. 
Other necessary measurements were the readings of the various heat meters, 
the temperature of the thermocouple reference plate, and wet and dry bulb 
temperatures. These measurements allowed determination of the local 
Stanton and Reynolds numbers. 

Transition from laminar to turbulent flow, which was stimulated by 
a strip of cellophane tape, occurred in each case on the first strip, as 
was indicated by the nature of the heat-transfer measurements. 


Data Reduction 

The heat leak from the back side of the plate and from the ends of 
each strip was evaluated from the heat-meter readings. The radiation 
heat transfer from the front of the plate was calculated from the temper- 
attire at the center of the strip, and the conduction between strips was 
calculated from the> center temperatures and an estimated value of the 
thermal resistance between strips. A first approximation to the convec- 
tive heat transfer was then found by subtracting the various heat leaks 
from the measured power input. A first approximation to the heat-transfer 
coefficient was obtained by dividing the approximate convection heat 
transfer by the temperature difference at the center of the strip. This 
allowed estimation of the ratio of the mean temperature difference to the 
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center temperature difference from the analysis of transverse tempera- 
ture distribution. Using the estimated mean temperature difference, sec- 
ond approximations to the radiation (based on the mean temperature), the 
heat transfer, and the conductance were calculated; then a new estimate 
of the mean temperature correction was made, again using the analysis. 
This iterative process was repeated until the heat- transfer rate, the 
mean temperature difference, and the heat- transfer coefficient were all 
in agreement. Usually only one iteration was required. 

Reduction of dynamic-pressure measurements in the usual manner led 
to calculation of the mass velocity. A humidity correction was made in 
determining the air density. The local Stanton and Reynolds numbers were 
then calculated, basing all fluid properties on the free-stream tempera- 
ture. The temperature -dependent fluid-properties correction, (T^/T^)®* 4 , 
was evaluated with the mean temperature of the strip, and the corrections 
on the Stanton numbers were of the order of 2 percent. 


Pressure Distribution and Free-Stream Turbulence Level 

After installing the plate in the wind tunnel, the angle of attack 
of the plate was adjusted to give the best possible pressure distribu- 
tion on the heated surface. The final pressure distribution (fig. 6) is 
considered satisfactory. 

The free-stream turbulence level was measured by hot-wire-anemometry 
techniques. The equipment used in these measurements is described in 
reference 15, and the results are shown by figure 7. The turbulence 
level was higher near the front of the plate and increased with increas- 
ing velocity. The turbulence intensity over most of the plate during 
typical runs was of the order of 2 to 3 percent. 


Precision 

An analysis of the experimental uncertainty indicated that the prob- 
able error in the local Stanton numbers is +3 percent, and the probable 
error in the local Reynolds number is ±1 percent. These values agree 
with the standard deviation of the experimental Stanton numbers from the 
isothermal equation (9), which was calculated as 4.5 percent. The anal- 
ysis of experimental uncertainty is presented in detail in appendix A of 
reference 1. 
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RESULTS AND DISCUSSION 
Heat-Transfer Data 

Eight isothermal heat-transfer runs were made at velocities ranging 
from 43 to 127 feet per second. Figure 8 shows the results of these 
tests, the data of which are tabulated in table I (a). The measured local 
Stanton numbers, corrected for fluid-property variations by the factor 

CT./’U are plotted against the local Reynolds numbers , ‘ based on the 

free- stream temperature. The data are compared with the von Kalman anal- 
ogy, employing the Schultz -Grunow friction factor (eq. (6)) and the sim- 
ple power equation (9) . The data are in good agreement with these rela- 
tions over the Reynolds number range 10^ < Re x < 3.5XL0®* The standard 
deviation of the experimental Stanton numbers from the power formula is 
4.5 percent. 


Velocity Survey Data 

Velocity surveys were taken at three points on the plate and at two 
different free-stream velocities. The purpose of these surveys was to 
demonstrate that the profiles were in agreement with the standard 
boundary-layer profiles and that thus there were no flow anomalies. 

These data are compared with the universal velocity profile for a flat 
plate (eqs. (15) to (18)) in figure 1. In reducing the survey data to 
a u + against y + basis, the wall shear stress was computed from the 
Schultz -Grunow expression for the friction factor (eq. (5) ) . The agree- 
ment of these data with the universal profile is quite good, except in 
the turbulent wake region where the correlation is not expected to be as 
satisfactory. The data are compared with the universal velocity defi- 
ciency profile of Schultz -Grunow in figure 9. Since the Schultz -Grunow 
friction-factor relation (5) was obtained from this deficiency profile, 
and since the data based on wall shear values from equation (5) agree 
quite well with the deficiency profile, it is felt that equation (5) is 
entirely adequate for determining the friction factors for the apparatus 
of the present investigation. 

The data are compared with the power profiles in figure 10, which 
shows that, in this Reynolds number range, the l/5. 6-power profile 



is the best fit of the data. The l/7 -power profile does not agree well 
in this Reynolds number range. The Schultz -Grunow velocity profile for 
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a Reynolds number of 10 is in good agreement with the data and with the 
l/ 5 . 6-power profile . The velocity survey data are tabulated in table 
1(b) . 

In reducing the velocity survey data, the hydrodynamic boundary- 
layer thickness 8 was determined by the method described earlier. The 
dimensionless velocity u/u was plotted against the distance from the 

wall on log-log paper, and the best "power formula fit" was obtained. 

The profile was then integrated to obtain the displacement thickness 5*, 
which can be related to 8. The boundary-layer thickness determined in 
this manner is the same one that is used in integral analysis $ this 
thickness is approximately the same as the 99-percent thickness (the y 
for which u/u^ = 0.99). 


Temperature Survey Data 

Temperature surveys in the boundary layer were made at three points 
on the plate for two different free-stream velocities (table l(c)). A 
typical profile is shown in figure 11. The purpose of these surveys was 
to obtain profiles for constant wall temperature that would serve as a 
basis of comparison for the profiles taken with a step wall-temperature 
distribution (ref. 2). As a matter of general interest, these data are 
also compared with the universal temperature profile predicted from the 
universal velocity profile. Because of the uncertainty in location of 
the thermocouple probe, it was necessary to "shift" the data so that the 
temperatures extrapolated to the wall value with the correct slope. The 
temperature gradient at the wall was determined from equation (9) and the 
Fourier heat-conduction equation 




-k 


~w 


which applies, since the flow near the wall is laminar. Both positive 
and negative shifts were necessary, indicating that the error was of a 
random nature. The method of the shift is indicated by figure 11. In 
no case did the shift exceed 0.002 inch, the average shift being about 
0.001 inch. 


The temperature survey data are compared with the universal temper- 
ature profile predicted earlier from the velocity profile in figure 2. 

In reducing the survey data to a t + against y + form, the wall shear 
stress was determined from the Schultz-Grunow friction factor (eq. (5)), 
and the wall heat flux was determined from the power relation (eq. (9)). 
The data agree well with the predicted profile in the laminar sublayer 
and in the buffer layer, but in the turbulent core the agreement is not 
too good. This is probably due to the fact that the eddy diffusivities 
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for heat and momentum are not equal in the turbulent core, as was as- 
sumed in the profile prediction. In fact, it is expected that the larg- 
est departure from this assumption occurs in the core. Reference 16 pro- 
poses an analysis from which the ratio of eddy diffusivities may be de- 
termined as a function of the Prandtl number and the quantity %/v. This 
analysis indicates that the ratio %/% is unity when Pr = 1 and 

£ m/ v = 1 > i>ht that e H /e M > 1 when e M /v > 1. Moreover, the departure 

from equal diffusivities is greatest when the turbulent portion of the 
momentum transfer is greatest, and this point occurs somewhere in the 
turbulent core. For air, the ratio of diffusivities is about 1.1 to 1.3. 
This explains why the data do not agree well with the predicted temper- 
atures in this region. 

Although the temperature and velocity profiles are not too "similar" 
when considered on a t + and u + against y + basis, the profiles do 
have similar shapes when the distances from the wall are characterized 
by the thicknesses of the respective hydrodynamic and thermal boundary 
layers. The velocity profiles were found to follow the l/5. 6-power 
formula. 



The temperature profiles can be represented by a similar relation, where 
the thermal boundary-layer thickness is employed: 



The temperature survey data are compared with power formulas in figure 
12. The l/5. 6-power curve is the best fit, indicating that the velocity 
and temperature profiles are indeed similar when each is referred to its 
respective boundary-layer thickness. This empirical observation is of 
considerable importance in integral analyses of turbulent-boundary-layer 
heat transfer. 

In reducing the temperature survey data, the thermal boundary-layer 
thickness bj was determined by the same method employed in evaluating 

the hydrodynamic boundary-layer thickness 6. The dimensionless temper- 
ature Q was plotted against the distance from the wall on log-log 
paper, and the "best power" behavior was determined. In each case it was 
found that 6 varied almost as the^di stance from the wall to the l/5.6 

power. The conduction thickness 5rp was then determined by numerical 

* 

integration of the temperature profile; bj, is related to for power 

profiles by 
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Srp = (l + m)s£ (31) 

l/m 

if 0 = (y/S^i) ' . This technique allowed determination of the same Sy 

that is used in integral analyses; the thicknesses determined in this 
manner were approximately the same as. the 99-percent thicknesses (the y 
for which 0 = 0.99). 


SUMMARY OF RESULTS 


The heat-transfer data for an isothermal plate are in good agree- 
ment with the best available analyses over the Reynolds number range 

C C 

10 < Re x < 3.5x10 . For air, a satisfactory power representation of the 
more complicated von Karman analogy is 

^ .-0.4 


StPr 0-4 = 0.0296 Re“ 0> 





(9) 


The Stanton, Prandtl, and Reynolds numbers in this relation are to be 
evaluated at the free-stream temperature. 


The velocity profiles agree well with the well-known profiles of 
Schultz -Grunow, and, in the Reynolds number range of the tests, the data 
may be represented by a power formula 


/ \l/5.6 

i-ffl 


(29) 


The temperature profiles are similar to the velocity profiles if 
each is characterized by its respective boundary-layer thickness. The 
temperature profile in the Reynolds number range of the tests can be 
represented by 



Stanford University, 

Stanford, Calif., October 22, 1957. 
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TABLE I. - EXPERIMENTAL DATA SUMMARY 


(a) He at -transfer data 


xl ° 3 St 


°F I BtU 
KhrHaq ft 


(hr)(sq ft)(-Op). 


St 

. ' 0. J 

*10 3 

25 

■p 

w 


xip3j 


127 ft/sec; t 


86.1° F; p o’o.C 

729 lb/cu ft I 

Nominal 

m 1 

20 ft/sec; t 

=■ 73.9° F; p 

a 0.0746 lb/cu 

483 


21.7 2.73 

2.77 

0.255 

32 .0 

22.8 

467 

20.5 

2.67 

2.72 

442 


19.4 2.41 

2.45 

.423 

32.1 

23.1 

431 

18.7 

2.42 

2.46 

386 


17.2 2.13 

2.17 

.580 

32.1 

22.9 

382 

16.7 

2.17 

2.20 

384 


17.1 2.11 

2.15 

.736 

32.1 

23.2 

381 

16.4 

2.14 

2.17 

-375 


16.7 2.06 

2.10 

.889 

32.1 

23.0 

369 

16.0 

2.08 

2.16 

367 


16.3 2.02 

2.06 

1 .045 

32.2 

23.3 

367 

15.8 

2.04 

2.08 

358 


15.8 1.97 

2.00 

1.196 

32.2 

23.3 

356 

15.3 

1.98 

2.01 

370 


16.2 2.01 

2.05 

1.353 

32.2 

23.3 

360 

15.5 

1.99 

2.03 

339 - 


14.7 1.85 

1.88 

1.507 

32.2 

23.3 

328 

14.0 

1.81 

1.84 

332 


14.3 1.79 

1.82 

1.661 

32.2 

23.6 

323 

13.7 

1.77 

1.80 

341 


14.8 1.84 

1.87 

1.823 

32.3 

23.5 

336 

14.3 

1.85 

1.88 

.326 


14.3 1.78 

1.82 

1.970 

32.2 

23.4 

323 

13.8 

1.78 

1.81 

.325 


14.3 1.79 

1.82 

2.13 

32.2 

23.4 

325 

.13.9 

1.80 

1.83 

317 


13.8 1.73 

1.76 

2.28 

32.1 

23.9 

321 

13.4 

1.74 

1.77 

320 


13.9 1.74 

1.77 

2.44 

32.2 

23.9 

326 

13.6 

1.76 

1.79 

325 


14.0 1.75 

1.78 

2.60 

32.3 

23.6 

320 

13.. 6 

1.75 

1.78 

317 


13.8 1.72 

1.75 

2 . 75 

32.2 

24.0 

315 

13.1 

1.70 

1.73 

314 


13.5 1.68 

1.71 

2.90 

32.1 

24.0 

316 

13.1 

1.69 

1.72 . 

322 


15.8 1.73 

1.76 

3.05 

32.1 

24.2 

318 

13.2 

1.71 

1.74 

308 


13.3 1.67 

1.70 

3.18 

32.1 

24.1 

303 

12.6 

1.64 

1.66 

289 


12.3 1.54 

1.59 

3.36 

32.1 

24.5 

287 

11. 7 

1.52 

1.55 

299 


12.7 1.58 

1.61 

3.51 

; 32.1 

24.5 

298 

12.2 ' 

1.58 

1.60 

113 ft/sec; 

fc « = 

74.7° F; p^ = 0.1 

0746 lb/cu 

ft 

Nominal i 

m l 

09 ft/sec; t 

« 74.9° F; p^ 

a 0.0 

745 lb/cu 


18.6 

18.0 

16.2 

15.6 
. 15.5 

14.8 

14.7 

14.9 
13.5 

14.3 
13.5 

13.4 
13.3 

2.56 

2.48 

2.23 

2.13 

2.12 

2.03 

2.01 

2.05 

1.84 
1.95 

1.85 
1.83 
1.81 

2.60 

2.53 

2.27 

2.17 

2.16 

2.07 

2.05 

2.09 

1.88 

1.99 

1.88 

1.86 

1.84 

0.235 
.386 
.528 
.676 
.816 
.960 
1.108 
1.247 
1.398 
1.539 
1.682 
1.828 
1 .978 

28.9 

29.2 
29.0 

29.2 

29.3 
• 29.3 

29.3 
•29.3 
29 ,4 
29 .4 

29.3 

29.4 
29.3 

23.2 
23.1 
.23.0 

23.5 

23.3 

23.6 

23.7 
23,5 

23.7 
23.9 
24.0 

23.8 
23.8 

387 

405 

362 

359 

348 

347 

335 

341 

309 

306 

318 

306 

309 

12.9 

1.77 

1.80 

2.11 

29.3 

24.4 

303 

13.0 

1.78 

1.81 

2.25 

29.3 

24.4 

309 

12.7 

1.74 

1.77 

2.40 

29.1 

24.3 

305 

12.5 

1.72 

1.75 

2.54 

29.3 

24.6 

296 

12.4 

1.71 

1.74 

2.68 

29.3 

24.6 

291 

12.5 

1.72 

1.75 

2.80 

29.0 

24.7 

298 

12.0 

1.66 

1.69 

2.95 

29.0 

24.5 

285 

11.1 

1.53 

1,56 

3.10 

29.1 

24.7 

267 

11.5 

1.58 

1.61 

3,25 

29.1 

24.8 

275 

66.3° F; 

p^ = 0.0757 lb/cu ft 

Nominal 

U = 

83 ft/s« 

18.5 

2.87 

2.92 

• 0.211 

22.6 

20.0 

328 

16.2 

2.51 

2.55 

.347 

22.6 

20.3 

284 

14.5 

2.21 

2.25 

.475 

22.6 

20.4 

248 

13.9 

2.15 

2.19 

.603 

22.7 

20.4 

255 

13.6 

2.11 

2.14 

.733 

22,6 

20.5 

249 

13.6 

2.10 

2.13 

.862 

22.7 

20.5 , 

243 

12.9 

1.99 

2.02 

.988 

22.7 

20,5 

233 

13.2 

2.05 

2.09 

1.115 

22.7 

20.5 1 

237 

12.1 

1.87 

1.90 

1.250 

22-7 

20.7 

224 

11.7 

1.81 

1.84 

1.373 

22.7 

21.1 

220 

12.3 

1.91 

1.94 

1 .502 

22.7 

20.9 

228 

11.6 

1.80 

1.83 

1.636 

22.7 

21.2 

223 

11.8 

1.83 

1-86 

1.729 

22.7 

21.0 

218 

11.4 

1.77 

1.79 

1.888 

22.7 

21.2 

214 

11.4 

1.76 

1.79 

2.02 

22.7 

21.2 

214 

11.3 

1.75 

1.78 

2.15 

22.7 

21.0 

210 

11.1 

1.72 

1.75 

2.27 

22.7 

21.1 

206 

10.9 

1.69 

1.73 

2.40 

22.7 

21.1 

208 

11 .3 

1.75 

1.78 

2.52 

22.6 

21.0 

205 

10.6 

1.66 

1.68 

2.65 

22.6 

21.1 

202 

9.92 

1 . 55 

1.58 

2.79 

22,6 

21.5 

190 

10.3 

1.60 

1.63 

2.92 

22.7 

21.5 

195 

70.6° F; 

a 0.0759 lb/cu ft 

Nominal 

U = 

43 ft/s€ 


■ 0.0760 lb/cu ft 



p .* 0.0765 lb/cu ft 





TABLE I. - Concluded. EXPERIMENTAL DATA SUMMARY 

(b) Velocity survey data 


Re x = 0.970*10 6 


Re x » 1 .210*10® 

Re x - 1, 738*10® 

Re x = 2.32*10 6 

= 40.9 ft/sec 


u a = 40.4 ft/sec 

« 99.6 ft/see 

u < = 99.1 ft/sec 

p” = 0.0771 lb/ cu 

ft 

p” = 0.0782 lb/cu ft 

p” » 0.0755 lb/cu ft 

p“ = 0.0757 lb/cu ft 


Velocity ratio, u/u 


(c) Temperature survey data 

Re -=0.729 xlQ 6 : .|Re„ .« 0.978xl0 6 |Re Y = 1.178xl0 6 |rc = 1.72*10 6 |Re T « 2.31xl0 6 ]Re = 2.78xl0 6 ’} 


t w - t = 27.3° F 


T> 

in. 

9 

in. 

e 

yv 

in. 

e 

T» 

in. 

e 


Q 

0 

0 

0 

0 

0 

0 

.004 

.128 

.004 

.137 

.005 

,129 

.003 

.242 

.005 

.178 

.005 

.181 

.006 

.177 

,004 

.278 

.006 

.219 

.006 

.225 

.007 

.217 

.006 

.363 

.007 

,252 

.007 

.267 

.008 

.245 

.008 

.443 

.008 

.277 

.008 

.306 

,009 

.282 

.010 

.476 

.009 

.314 

.009 

.326 

.010 

.306 

.012 

.506 

.010 

.347 

.010 

.358 

.011 

.325 

.014 

. 524 

.011 

.364 

.011 

.347 

.012 

.346 

.016 

,535 

.013 

.409 

.013 

.410 

.014 

.370 

.018 

.549 

.015 

.442 

.015 

.438 

.016 

.420 

.023 

.575 

.017 

.459 

.017 

.463 

.018 

.442 

.028 

.590 

.019 

.492 

.019 

.486 

.020 

.462 

.033 

.608 

.021 

.512 

.021 

,506 

.024 

.495 

.038 

,619 

.023 

.521 

.023 

.524 

.029 

.523 

.043 

.630 

.028 

.550 

.028 

.544 

.034 

.551 

,048 

.642 

.033 

.575 

.033 

.577 

.044 S 

.579 

.058 

.664 

.043 

.604 

.043 

.600 

.054 

,604 

.068 

.678 

.053 I 

.633 

.053 

.612 

.064 | 

.615 'I 

.078 

.689 

.063 

.649 

.063 

.631 

.074 

.631 | 

.098 

.714 

.073 

.666 

.073 

.651 

.084 

.648 

.148 

, 754 

.083 

.678 

.063 

,660 

.104 

.668 

,198 

.795 

.103 1 

.703 

.103 

.680 

.154 

.704 

.248 

.827 

.153 1 

.748 

.153 

,720 

.204 

. 728 

.298 

.865 

.203 ! 

.790 

.203 

.775 

.254 

. 760 

.348 

.890 

.253 

.823 
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(b) Plan view showing heated- section dimensions (in inches). 
Figure 3. - Concluded. Experimental flat plate. 
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Concluded. Wind-tunnel installation 
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Figure 5. - Thermocouple probe. 
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